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Abstract 

We investigate the presence of irrelevant operators in the 2d Ising 
model perturbed by a magnetic field, by studying the corrections in- 
duced by these operators in the spin-spin correlator of the model. To 
this end we perform a set of high precision simulations for the corre- 
lator both along the axes and along the diagonal of the lattice. By 
comparing the numerical results with the predictions of a perturbative 
expansion around the critical point we find unambiguous evidences of 
the presence of such irrelevant operators. It turns out that among 
the irrelevant operators the one which gives the largest correction is 
the spin 4 operator -|- which accounts for the breaking of the 
rotational invariance due to the lattice. This result agrees with what 
was already known for the correlator evaluated exactly at the critical 
point and also with recent results obtained in the case of the thermal 
perturbation of the model. 
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1 Introduction 



Since the seminal work of Belavin, Polyakov and Zamolodchikov [0] much 
progress has been done in the description of 2d statistical models at the 
critical point. In particular in all the cases in which the critical point theory 
can be identified with one of the so called "minimal models" a complete list of 
all the operators of the theories can be constructed. Moreover, it is possible 
to write differential equations for the correlators and in some cases find their 
explicit expression in terms of special functions. Following these remarkable 
results, in these last years a lot of efforts have been devoted to extend them 
also outside the critical point. This can be achieved by perturbing the CFT 
with one (or more) of its relevant operators. Much less is known in this case. 
The most important result in this context is that for some particular choices 
of CFT's and relevant operators these perturbations give rise to integrable 
models @, ^. In these cases again we have a rather precise description of the 
theory. In particular it is possible to obtain the exact asymptotic expression 
for the large distance behavior of the correlators 0]. From this information 
several important results (and in particular all the universal amplitude ratios) 
can be obtained. 

Similar efforts have also been devoted in trying to study vacuum expecta- 
tion values, amplitude ratios and correlators involving irrelevant operators. 
However progress in this direction is limited by the lack of results (both nu- 
merical and analytical) from actual statistical mechanics realizations of the 
models. 

The only notable exception to this state of art is the 2d Ising model. 
In this case, thanks to a set of remarkable results on the 2 and 4 point 
correlators §|, 01 several interesting results on the contribution of irrelevant 
operators have been obtained both in the case of the model at its critical 
point 1^ and for T 7^ Tc (i.e. for the theory perturbed by the energy 

operator). Very recently these results have been used to study up to very high 
orders the contribution of irrelevant operators to the magnetic susceptibility 
for T ^ Tc fry]. In this paper we shall try to extend this analysis to the 



case of the magnetic perturbation of the Ising model. In particular, we 
shall look at the contributions due to irrelevant operators to the spin-spin 
correlator in presence of an external magnetic field. In particular we shall 
mainly concentrate on the energy momentum tensor TT, and the operator 
_l_ 2^2 appear as a consequence of the breaking due to the lattice of 
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the full rotational symmetry) which are the most important (i.e. those with 
smallest (in modulus) renormalization group eigenvalue) among the irrelevant 
operators. 

In this case there is no exact expression for the correlators or for the free 
energy of the model and the only way that we have to judge of the validity 
of the CFT predictions is to compare them with the results of numerical 
simulations. Notwithstanding this there is a number of good reasons to 
choose exactly this model and this observable. Let us look at them in detail. 

• The Ising model in a magnetic field is known to be exactly integrable 
thus, even if it is not exactly solvable, a lot of important informations 
can all the same be obtained, in particular very precise large distance 
expansions exist for the correlators and the vacuum expectation values 
of the relevant operators are known in the continuum limit [|ll|]. 

• The model is an optimal choice from the point of view of numerical 
simulations, since very fast and efficient algorithms exist to study it. 

• The model can also be realized as a particular case of the so called 
IRF (Interaction Round a Face) models . In this framework several 
interesting results can be obtained on the spectrum of the model |]13| . 

• By looking at the spin-spin correlator instead of the susceptibility we 
may study the effects of the non zero spin operators which appear as 
a consequence of the breaking of the rotational symmetry due to the 
lattice (see below for a precise definition). These operators also appear 
in rotational invariant quantities like the susceptibility, but only at the 
second order i.e. at such an high power of the perturbing 
constant that they can be observed only if the exact solution of the 
model is known (as in the papers 0, [1^]) and cannot be seen if one can 
only resort to numerical simulations. 

• A very interesting result of ||10[ (which had already been anticipated 
twenty years ago by Aharony and Fisher |TB[ and also confirmed in ) 



is that in the thermal perturbation of the Ising model the energy mo- 
mentum tensor TT seems to be absent. The first correction which 
involves irrelevant operators is thus due to the spin 4 operator + T^. 
This result should also hold at the critical point where it can be inde- 



pendently observed by using transfer matrix methods [0. However, 



2 



apparently, it disagrees with the explicit form of the correlator at the 
critical point (which is exactly known on the lattice) which contains a 
scalar correction which can only be interpreted as due to the TT term. 
One of the aims of the present paper is to understand this puzzle. 

If one is interested in studying irrelevant operators it is mandatory to 
look at the short distance behaviour of the correlator. In this respect 
the natural framework in which one must operate is the so called IRS 
expansion |T8|, |19|. This approach has been recently discussed in great 
detail exactly in the case of the Ising model in a magnetic field in 



which we are interested here. It is only by using the results of ||2^ as 
input of our analysis that we shall be able to reach the high level of 
precision which is needed in order to observe the very small corrections 
which are the signatures of the irrelevant operators. 

This paper is organized as follows. In sect. 2 we shall briefly summarize 
some known results on the 2d Ising model in a magnetic field. In sect. 3 we 
shall discuss the most important contributions due to the irrelevant operators 
to the spin-spin correlators and shall evaluate their magnitude and behaviour. 
In sect. 4 we shall present the numerical simulations that we have performed 
and finally in sect. 5 we shall discuss the comparison between numerical results 
and theoretical predictions. Sect. 6 is devoted to some concluding remarks. 



2 The 2d Ising model in a magnetic field. 

We shall be interested in the following in the Ising model defined on a 2d 
square lattice of size L with periodic boundary conditions, in presence of an 
external magnetic field H. The model is defined by the following partition 
function: 

Z = J2 e'^^^i''^)"'"''^"^''"''^ (1) 

where the notation denotes nearest neighbour sites in the lattice. In 
order to select only the magnetic perturbation, j3 must be fixed to its critical 
value: 

/3 = = - log (V2 + 1) = 0.4406868... 
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by defining hi = (3cH we have 

The magnetization M[h) is defined as usual: 

^W-^^(log Z)\,=,^ = {^Y.<^,). (3) 

where N = L"^ denotes the number of sites of the lattice. 

Eq. (ID is the typical partition function of a perturbed critical model. With 
the choice P = Pc the only perturbing operator is 

^^ = ^E^* ' (4) 

i 

We shall call in the following ai as the spin operator (more precisely the 
lattice discretization of the spin operator). Notice that the mean value of ai 
coincides with M{h): 

{ai) = M{h) (5) 

Our goal in this paper is to study the contribution of the irrelevant opera- 
tors to the spin-spin correlator. To this end we shall first study the model 
at the critical point (sect. p.l|) , we shall then switch on the magnetic field 
(see sect.|2.2|) and discuss the modifications that it induces in the spin-spin 
correlator 

2.1 The Ising model at the critical point 
2.1.1 Operator content. 

The Ising model at the critical point is described by the unitary minimal 
CFT with central charge c = 1/2 [|l|. Its spectrum can be divided into 
three conformal families characterized by different transformation properties 
under the dual and Z2 symmetries of the model. They are the identity, spin 
and energy families and are commonly denoted as [/], [a], [e]. Each family 
contains a "primary" field (which gives the name to the entire family) and an 
infinite tower of "secondary" field (see below). The conformal weights of the 
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primary operators are hj = 0, hfj = 1/16 and h^ = 1/2 respectively. Thus we 
see that in the Ising model the set of primary fields coincides with that of the 
relevant operators of the spectrum (remember that the relationship between 
conformal weights and renormalization group eigenvalues is: y = 2 — 2h). 
This is a peculiar feature of the Ising model only, and is not shared by any 
other minimal unitary model. Thus in this case the irrelevant operators 
are bound to be secondary fields. Since in this paper we are particularly 
interested in the irrelevant operators, let us study in more detail the structure 
of the three conformal families. 

• Secondary fields 

All the secondary fields are generated from the primary ones by apply- 
ing the generators L_j and L_j of the Virasoro algebra. In the following 
we shall denote the most general irrelevant operator in the [a] family 
(which are odd with respect to the Z2 symmetry) with the notation 
and the most general operators belonging to the energy [e] or to the 
identity [/] families (which are Z2 even) with and rji respectively. It 
can be shown that, by applying a generator of index k: L_k or L_k to 
a field (f) (where (t) — I,e,a depending on the case), of conformal weight 
^0, a new operator of weight h = + k is obtained. In general any 
combination of L_j and L_j generators is allowed, and the conformal 
weight of the resulting operator will be shifted by the sum of the indices 
of the generators used to create it. If we denote with n the sum of the 
indices of the generators of type L_j and with n the sum of those of type 
L_i the conformal weight of the resulting operator will be /i^ + n + n. 
The corresponding RG eigenvalue will he y — 2 — 2h^ — n — n. 

• Nonzero spin states 

The secondary fields may have a non zero spin, which is given by the 
difference n — n. In general one is interested in scalar quantities and 
hence in the subset of those irrelevant operators which have n = n. 
However on a square lattice the rotation group is broken to the fi- 
nite subgroup C4 (cyclic group of order four). Accordingly, only spin 
0, 1, 2, 3 are allowed on the lattice. If an operator of the continuum 
theory has spin j G N, then its lattice discretization 0; behaves as a 
spin j (mod 4) operator with respect to the C4 subgroup. As a conse- 
quence all the operators which in the continuum limit have spin j = 4N 
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with non-negative integer can appear in the lattice discretization of 
a scalar operator. 

• Null vectors 

Some of the secondary fields disappear from the spectrum due to the 
null vector conditions. This happens in particular for one of the two 
states at level 2 in the a and e families and for the unique state at 
level 1 in the identity family. From each null state one can generate, 
by applying the Virasoro operators a whole family of null states hence 
at level 2 in the identity family there is only one surviving secondary 
field, which can be identified with the stress energy tensor T (or T). 

• Secondary fields generated by L i 

Among all the secondary fields a particular role is played by those gen- 
erated by the L_i Virasoro generator. L_i is the generator of trans- 
lations on the lattice and as a consequence it has zero eigenvalue on 
translational invariant observables. 

• Quasiprimary fields 

A quasiprimary field \Q) is a secondary field which is not a null vector 
(or a descendent of a null vector) and satisfies the equation 

^iig) = o (6) 

These fields play a central role in our analysis since they are the only 
possible candidates to be irrelevant operators of the model. 

By imposing eq.(^ it is easy to construct the first few quasiprimary 
operators for each conformal family. For our analysis however the two 
lowest ones are enough. They both belong to the conformal family of 
the identity. Their expression in terms of Virasoro generators is: 

Q\ = ^-2\V) (7) 
Q\ = {^\-\i^^m (8) 

(we use the notation Q\ to denote the quasiprimary state at level n in 
the 1] family). 

From these fields we can construct two irrelevant operators, which both 
have conformal weight 4 and RG eigenvalue —2. 
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1] The first combination is QlQl which has spin zero and can be 
identified with the energy momentum tensor TT. 

2] The second combination is Q\ + Ql which has spin 4 and can be 
identified with the combination + T^. 

This last contribution appears as a consequence of the breaking of the 
full rotational invariance due to the lattice^. 



2.1.2 Structure constants. 

Once the operator content is known, the only remaining information which is 
needed to completely identify the theory are the OPE constants. The OPE 
algebra is defined as 

<f.W$,(o) = E4'^W*w(o) (9) 
{k} 

where with the notation {k} we mean that the sum runs over all the fields 
of the conformal family [k]. The structure functions C^j{r) are c-number 
functions of r which must be single valued in order to take into account 
locality. In the large r limit they decay with a power like behaviour 

Cf.(r) ~ |r|-'''™(^S) (10) 



whose amplitude is given by 



= Jim Cj.(r) \r\'-(^^.\ (11) 

Several of these structure constants are zero for symmetry reasons. These 
constraints are encoded in the so called "fusion rule algebra" which, in the 
case of the Ising model is. 

m = [1] 

[a][e] = [a] (12) 
MM = [1] + 



^Notice that if we would be interested in scalar quantities, this term would disappear 
even on the lattice at the first order and could contribute only at the second order. This is 
the case for instance of the susceptibility recently discussed in , in which this operator 
gives a contribution only at order and not t^. However in the present case since we are 
interested in a correlator, which defines a preferential direction on the lattice this term 



can contribute already at the first order. This fact will be discussed in detail sect. 2.1.5 
below. 
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By looking at the fusion rule algebra we can see by inspection which are the 
non- vanishing structure constants. 

The actual value of these constants depends on the normalization of the 
fields, which can be chosen by fixing the long distance behaviour of, for 
instance, the aa and ee correlators. In this paper we follow the commonly 
adopted convention which is: 

((T(x)cr(O)) = , |x| — i> CX) (13) 

(e(x)e(O)) = ^ , |x| oo. (14) 

fI 

With these conventions we have, for the structure constants among primary 
fields 

CI. = CI = CI, = (15) 
C'.V = C:,i = Cl = = 1 (16) 

and 

C,. = = \- (17) 
2.1.3 Continuum versus lattice operators. 

Our main interest in this paper is the spin-spin correlator on the lattice. This 
rises the question of the relationship between the lattice and the continuum 
definitions of the operator a. In the following we shall denote the lattice dis- 
cretization of the operators with the index /. Thus a denotes the continuum 
operator and ai the lattice one. 

In general the lattice operator is the most general combination of contin- 
uum operators compatible with the symmetry of the lattice one. If we are 
exactly at the critical point this greatly simplifies the analysis, since only 
operators belonging to the [a] family are allowed. Moreover (due to the pe- 
culiar null state structure of the spin family) the first quasiprimary operator 
in the spin family appears at a rather high level and can be neglected in 
our analysis. Thus as far as we are interested only at the critical point the 
spin operators on the continuum and on the lattice are simply related by a 
normalization constant which we shall call in the following R^j. The simplest 
way to obtain R^j is to look at the analogous of eq.(|13[) on the lattice. This 
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is a well known result which we report here for completeness. The large 
distance behaviour of the correlator at the critical point is 



Rl 



11/4 



(18) 



where rij denotes the distance on the lattice between the sites i and j and 

Rl = e3?'(-i)25/24 = 0.70338... (19) 
By comparing this result with eq.(|l3|) we find 

ai = = 0.83868.. .a (20) 



2.1.4 The lattice Hamiltonian at the critical point 

The last step in order to relate the continuum and lattice theories at the 
critical point is the construction of the lattice Hamiltonian (let us call it 
Hiat) at the critical point. As above, the lattice Hamiltonian will contain all 
the operators compatible with the symmetries of the continuum one. In this 
case all the operators belonging to the [a] family are excluded due to the 
Z2 symmetry. Also the operators belonging to the [e] family are excluded 
for a more subtle reason. The Ising model (both on the lattice and in the 
continuum) is invariant under duality transformations while the operators 
belonging to the [e] family change sign under duality, thus they also cannot 
appear in Hiat{t = 0). Thus we expect 

Hiat = HcFT + w° / d^^Vi, Til e [/] , (21) 

where He ft is the continuum Hamiltonian and the are constants. As 
mentioned above we can keep in this expansion only the first two terms 
which are respectively TT and T"^ + T'^. 



2.1.5 The spin-spin correlator on the lattice at the critical point. 



We have already presented above, in eq. (|Tq) the large distance expansion of 
the spin-spin correlator on the lattice at the critical point. 
Thanks to the exact results obtained in P, |^, 0] 



see 



22 



for a review) we 

have much more informations on this correlator at the critical point. In order 
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to discuss these results let us first introduce a more explicit notation for the 
spins. We shall denote in the following with aM,N the spin located in the site 
with lattice coordinates (M, A^). By using the translational invariance of the 
correlator we can always fix one of the two spin in the origin. Thus the most 
general correlator can be written as: (o"o,oC"M,Ar)- A particular role will be 
played in the following by the correlator along the diagonal: (cro,ocrAr,Ar) and 
the one along the axis (o"o,oC"o,Af)- 

Among the various results on the critical correlators two are of particular 
relevance for us: 

• Remarkably enough, exact expressions exist for the correlator both 
along the diagonal and along the axis for any value of N . These ex- 
pressions are rather cumbersome and we shall not report them here. 



They can be found in p2|. 



• An asymptotic expansion exists for large values of the separation rij. 
For a square lattice, the first three terms have been explicitly evaluated 
ini. 

This expansion takes the following form (for further details see |]^) 
log((To,ocrM,iv) = log A - ^ log r + Ai (6')r + 

-4 



+ A2(%-' + A3(%-^ + 0(r-«) (22) 



where 



A^{e) = 2-S(-l + 3cos40) 

A2{9) = 2"^^(5 + 36cos4^ + 36cos80) 

As^e) = 3"^2"^^(-524- 324 cos 4^ + 24732 cos 80 + 28884 cos 12^X23) 

and 

= 1 (M2 + N^) 
M = r a/2 sin 6* 

N = r \/2cos0 . (24) 
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In the following we shall only be interested to the first correction since higher 
order corrections are beyond the resolution of our data. Hence we end up 
with 

{cro,ocrM,N) = ^ + 2"'(-l + 3 cos 4^) + . . .) (25) 

Remarkably enough this expansion perfectly agrees with what one finds by 
assuming the presence in the Hamiltonian of the model of the two irrele- 
vant operators discussed in sect. |2.1.I| . In fact both TT and + have 
RG eigenvalue -2, thus, if they are present, they should contribute to the 
correlator exactly with a term proportional to Moreover we expect 

that the + operator, which has spin 4, should give a term proportional 
to cos(46') while the scalar operator TT should give a contribution without 
9 dependence. This is exactly the pattern that we find in eq.(p5|) (see for 
instance pag. 218 of for a detailed discussion of this point). However 
this remarkable agreement also rises a non trivial problem. In fact the high 
precision analysis of |T^, fl^ clearly exclude the presence in the lattice 



Hamiltonian of the TT operator (while they both confirm that the + 
is indeed present) It is thus not clear which could be the origin of the scalar 
term in eq.(^). A possible solution to this puzzle is to notice that the results 
of eq. (|25|) and those of of [10|, |14| are obtained with two different choices 



of the coordinates of the 2d plane. 

In fact eq.(p5D is written in terms of the "continuous" variable r (which, is 
the one that we must choose if we want to match the lattice results with the 
continuum limit ones). Both the results of and [|14| are instead obtained 



in the "lattice reference frame" (which is the most natural variable on the 
lattice) in which there is no a/2 when comparing the distances along the axes 
and along the diagonals. Hence, to compare the correction to the spin-spin 
correlation function with the findings of of [|I0|, [|I^ we must rewrite (p5|) 
it in terms of lattice coordinates. This can be easily performed by using the 
relations (24). 



After some algebra we obtain 
• axis correlator 

A 2i/« / 1 



acoao.iv) = [1 + -^ + OiN-') ) (26) 
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• diagonal correlator 

Wo,oaN,N) = ]^ (l - ^ + OiN-')) ■ (27) 

These same expansion can also be obtained by directly looking at the exact 
lattice results for the correlators (see for instance where these terms 
are obtained in full detail). 

Looking at eqs.(^, ^) we see that with the lattice choice of coordinates 
the term exactly changes its sign as we move from the axis to the 

diagonal. This is exactly what one would expect for the contribution of a 



spin 4 operator, and thus it is apparent from eqs.(p6|, |27D that no scalar 
correction appears at order l/N"^ in perfect agreement with the results of 

ofjig, 0. 



It is only the change of coordinates of eq.(p4D which induces in the con- 
tinuum limit a scalar term (which we may well identify in this limit with a 
TT type contribution). 

2.2 Adding the magnetic field. 
2.2.1 The continuum theory. 

The continuum theory in presence of an external magnetic field is represented 
by the action: 

A = Ao + h J d^x a{x) (28) 

where is the action of the conformal field theory. 
As a consequence of the applied magnetic field the structure functions 
acquire a h dependence so that we have in general 

$i(r)%(0)=Ecf (^,^)"^'W(0) • (29) 

{k} 

Also the mean values of the a and e operators acquire a dependence on 
h. Standard renormalization group arguments (for un updated and thor- 
ough review on renormalization group theory applied to critical phenomena 
2^1) allow one to relate this h dependence to the scaling dimensions of 



see 



the operators of the theory and lead to the following expressions: 

{a)h = A^hTs + ... (30) 
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{e)h = AM-^+... (31) 



The exact value of the two constants A^j and can be found in |25| and p6| 
respectively 

A„ = — . 5- — = 1.27758227.. , (32) 

15 sin ^ + sin ^ + sin ' ^ ^ 



with 



and 



4sin^rfi) /47r2rf|)r2(-^ 




r (!) r (^) V r (i) p (a) 

A, = 2.00314... . (34) 



2.2.2 Continuum versus lattice operators. 

In presence of a magnetic field also the fields belonging to the energy and 
identity families can appear in the relation between the lattice and the con- 
tinuum version of the spin operator. The most general expression is 

(Ti = foihi)a + hif^{hi)e + f[{hi)cr, + hifl{hi)e, + hif,\hi)r]„ i G N (35) 

where f[{hi) flihi) and /-{hi) are even functions of hi. With the notation 
(Tj, and rji we denote the secondary fields in the spin, energy and identity 
families respectively. This is the most general expression, however, as a 
matter of fact, only the first term is relevant for our purposes (all the higher 
terms give negligible contributions) 

ai = RrjG + hpie (36) 

The determination of pi is rather non-trivial. We shall discuss it in the 
appendix. It turns out that 

Pi ~ 0.0345 (37) 

It is also important to stress that the lattice and continuum values of the 
magnetic field do not coincide but are related by 

hi = Rhh (38) 

with Rh = 1.1923.. (see [BO] for details). 
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2.2.3 The lattice Hamiltonian. 



As it happened for the conversion from a to cfi also in the construction of 
the perturbed Hamiltonian on the lattice, several new operators belonging 
to the energy and spin family must now be taken into account. However 
it turns out that all the quasiprimary fields of these two families appear at 
rather high level and can be neglected in the present analysis. Thus, as far 
as we are concerned, the only effect of switching on the magnetic field is that 
the constants in front of the TT and + terms in eq.(|21|) acquire an h 
dependence. For symmetry reasons these must be even analytic functions of 
h. 



2.2.4 The spin-spin correlator. 

The short distance behaviour of the spin-spin correlator in presence of a 
magnetic field can be obtained by using the so called IRS approach. A 



detailed discussion of this approach can be found in [jT8|, [T^ the particular 
application to the Ising correlators are discussed in |^ to which we refer for 
details. We only list here the results. Setting t = \h\ |r|^^/^ we have up to 

o{e) 

(a(0)a(r))|r|V^ = Bl^ + BIj^/^^ + BIj^^'^' + 0{e) (39) 

with 

^ 1 





= ci„ = 


BL 


= A^Cl^ - 


bL 





1.00157... 
= -0.51581... 



(40) 



This result holds in the continuum theory. By using the known conversion 
between continuum and lattice units we can obtain the spin-spin correlator 
on the lattice: 

{a^iOMr))\n,\'/' = Bl, + Bl//'' + Bl.tf'' + 0(1^) (41) 

Where is the distance between the two spins on the lattice, measured 
in units of the lattice spacing, ti is defined as ti = \hi\ |rjjp/^ and 

Bl, = 0.703384... 
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Bl^^i = 0.641409... 
Bl^^ = -0.300749... 

(42) 

In the next section we shall see which contributions must be added to 
this expression as a consequence of the irrelevant fields. 

3 Contribution from irrelevant operators to 
the spin-spin correlator in presence of a 
magnetic field. 

There are two main sources of contributions due to irrelevant operator to 



eq.(|41|) above. They have different origin and contribute in a different way 
to the spin-spin correlator. The first correction can be understood as a per- 
turbative contribution to the spin-spin correlator due to irrelevant operators. 
We shall discuss it in subsec. 3.1. The second correction arises from the new 
terms in the relation between ai and a, due to the presence of a magnetic 
field. We shall discuss it in subsec. 3.2. 



3.1 Irrelevant operators in the perturbed Hamiltonian. 

The contributions to the scaling function due to the (possible) presence of 
additional irrelevant operators in the lattice Hamiltonian can be studied by 
means of standard perturbative expansions around the CFT fixed point. 
In order to address the problem in presence of an external magnetic field we 
must consider the CFT pertaining to the critical Ising model as perturbed by 
a mixed relevant /irrelevant perturbation. Hence there will appear correction 
terms to the scaling behaviour which both depend on the external magnetic 
field and show a non-trivial (i.e. non scaling) dependence on the spin-spin 
distance Vij. Moreover if the irrelevant operator in which we are interested 
breaks the rotational invariance (and this is the case for instance of the 
(T^ _l_ j-z-j QpQj-g^^QY discussed above) we must expect a dependence on the 
angle 6 between the principal axes of the lattice and the direction of the 
correlator. General symmetry arguments and dimensional analysis strongly 
constrain these terms. If we keep in our analysis only the first two irrelevant 
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operators TT and (T^ + T^) the most general expression for the scahng 
function turns out to be: 



'V • ■ \ \ 'Y' ■ ■ \ 



+ BL,fl + #^V;"" + 0(if) (43) 



where ai, 02 are unknown constants and b{6), c{9) are unknown functions of 
the form 

b{e) = h + b2 cos{4e), (44) 

c{e) = C1 + C2 cosine), (45) 

in the angular variable 6. 

Let us analyze the above corrections in detail. 

• Correction to Bl • 



(T(J, 



At the lowest order in the irrelevant coupling, we can identify a scalar 
correction proportional to and a spin 4 correction proportional 

to cos(4^) Both these terms acquire a h dependence due to the 

magnetic relevant perturbation. The most general expression, compat- 
ible with the symmetries of the model is: 



where Pi{h) and P2{h) are even polynomials in the magnetic field h 
which in the /i — > limit must agree with the asymptotic expansion 
reported in eg. (p5D . Expanding the polynomials up to the first term in 
h we find 

P2{h) = 3 + a2h\ (47) 



from which the first term in eq.(^) follows. It is easy to see that the 
terms proportional to Oi and 02 are highly suppressed due to the h'^ 
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power. They behave as the 0{t^) that we have systematically neglected 
in the previous section. As a matter of fact we have not been able to 
see these terms in our numerical data and we shall neglect them in the 
following. 

Corrections to B^^^ and Bl^f. 

Following the same line discussed above we find in this case: 



where b{6) and c{9) are the most general mixture of scalar and spin 
4 terms. They can be expanded in powers of cos(4^^). keeping only 
the first two orders in the expansion we end up with the expressions 
reported in eq.s (0) and ( ^5D where &i, 62; Ci and C2 are unknown con- 
stants. 

While the term proportional to c{6) cannot be detected within the 
precision of our data (we shall further discuss this point at the end 
of the next section), the magnitude of the one proportional to b{9) is 
definitely larger than our numerical uncertainties. Thus we expect to 
be able to observe such a correction. 



3.2 Irrelevant operators in ai. 



We have seen above (see eq.(P^)) that in presence of a magnetic field the 
relation between lattice and continuum spin operator becomes more com- 
plicated and a term proportional to he appears. Strictly speaking e is not 
an irrelevant operator, however this is only an accident, the following terms 
(those that we neglect) in the correspondence between ai and a would indeed 
be irrelevant operators. Moreover, exactly as it happens for the irrelevant 
operators, also the he term gives a subleading contribution to the spin-spin 
correlator. For these reasons we have included also this correction in this 
paper. 



Inserting eq. (pq ) in the spin-spin correlator we find 

{c7i{0)ai{r)) = Rl{a{0)a{r)) + 2hR^pi{a{0)e{r)) + hyMO)e{r)) (49) 
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Since we are interested in keeping only terms below O(t^) the last term in 
eg ■ (|49|) can be neglected. Using the known result (see for details) 

(a(0)e(r))|r|9/« = Bl^f/^^ + Bj^t + 52^^23/15 ^ o(p/i5) 

with 



we can rewrite eq.(H9p as 



bL = 


^aC^^ — 


0.63879... 


Bl = 


d^Cl, = 


3.29627... 


Bl = 




= - 1.82085 


1) as 







+ 2hR^pM-"\Blt'l''> + Blt + Blf'l'') (51) 

Only the first term in the second line of eq. (|5l|) gives a contribution below 
Oit^) thus we end up with 

{<ym<yi{r)M" = Rl{Bl + BlJi^' + Blj^'i^') 

+ 2R,p,t\r\-'^''Blt^'^' (52) 

We see that the new term can be considered as a correction, proportional 
to to the t^^/^^ term of eq.(|3): 

{^m^i{r))\n,\'/' = Bl, + Bl„/I'' 

+ {RlBl + 2RMr\-^^/'Bl)/Rf'h]'/'' + 0(1^ 

Inserting the values of the various constants and of pi we end up with 

(-KO)aKr))|r.//^ = 5^., + 5^/^ + (^L, + «3|rr^^/«)tr^ + 0(tn (54) 

with as = 0.0307... 

Unfortunately this new contribution is too small to give a reliable signa- 
ture within the precision of our data. In fact the amplitude of the correction 
03 |r I "22/8^16/15 distances greater than one lattice spacing is always below 
(or at most of the order of) 10~^ and cannot be disentangled from the other 
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sources of corrections within the precision of our data which ranges from 
10~^ to 10~^. The situation is very similar to that of the correction propor- 
tional to c{6) discussed in the previous section which in fact has the same t 
dependence. 

However it is interesting to observe that both these corrections are al- 
most in threshold to be observed and it is well possible that with the next 
generation of simulation they could be detected. For this reason we included 
also their analysis in the present paper. 



4 The simulation 

We studied the model with a set of Montecarlo simulations using a Swendsen- 
Wang type algorithm, modified so as to take into account the presence of an 
external magnetic field. We used the same program and simulation setting 



of ref. |]20| to which we refer for a detailed discussion of the performances of 
the algorithm and of finite size effects. We simulated the model for twelve 
different values of the magnetic field, ranging from hi = 4.4069 x 10"'^ to 
hi = 8.8138 X lO-l 

For all the values of hi that we simulated, we studied the correlator: 
(cr(O)cr(r)), for r = 1, ■ ■ ■ , 10 both in the diagonal direction and along the 
axes. This is an important feature of the present analysis, since it is exactly 
by comparing the correlator along the axes versus the diagonal one that we 
can extract the spin of the perturbing operators. 

For all the simulations the lattice size was chosen to be L = 200 (the 
analysis of shows that for all the values of hi that we studied this is 



enough to avoid finite size effects). Each two measures of the correlators were 
separated by five SW sweeps. For each single correlator we performed 4 x 10^ 
measures. The values of hi that we studied are reported in tab.^ where we 
also listed for completeness the corresponding values of the correlation length 
^. This quantity is very important , since it defines in a quantitative way 
what we mean as "short distance behaviour" . Short distance means "smaller 
than the correlation length". This is in fact the regime in which we expect 
that the IRS analysis should hold and in which we may hope to observe 
effects due to the irrelevant operators. ^ is given (in lattice units) by: 

Cihi) = 0.24935.. ./i;"^ (55) 
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(see for details on this equation and on the continuum to lattice conver- 
sion of ^) . 



Table 1: The values of hi (measured in units of f3c) that we studied and the 
corresponding correlations lengths. 



hi 




0.001 


15.4 


0.002 


10.6 


0.003 


8.5 


0.004 


7.3 


0.005 


6.5 


0.006 


5.9 


0.007 


5.4 


0.008 


5.1 


0.009 


4.8 


0.01 


4.5 


0.015 


3.6 


0.02 


3.1 



5 Comparison between predictions and nu- 
merical results. 

Let us start this section by writing the expression that we expect for the 
spin-spin correlator in a magnetic field, according to the analysis performed 
in the sections 2 and 3 above. We have: 



{ai{0)aiir))\ri,\'/' = {aiiO)ai{r)),=o\ri,\'/' + (a, + ^ + . . .) t', 



where {ai{0)ai{r))h=o\rij\'^^^ is given by eq.(||); ai = Bj^i and bi = B^^i, 
while 02(6') and 62 (^) are in principle unknown functions of the angular vari- 



+ 




(56) 



able e. 
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In this section we shall compare this expression with our numerical estimates 
for the correlators along the diagonal and the axis. This will allow us to ob- 
tain a rather precise estimate of the function 02(6') in the two cases ^ = 
and 9 = 7r/4. With these two informations we shall also be able to study the 
spin content of the operators which generate this contribution. 
In the following we shall confine ourselves to the study of 02(6') since precision 
of our data is not enough to extend our analysis also to the function 62 (6*) • 
We performed our analysis in three steps. 

1] As a first step we constructed the combinations 

G^{r, h) = {ai{0)aiir)) - {ai{0)aiir))^=o (57) 

using the fact that the critical point correlator is exactly known. Thus 
we have: 

+ (6, + M^ + ...)t-/- + 0(t?). (58) 

2] Then we study, at fixed value of r, the hi dependence of G^{r, hi). Ac- 
cording to eq.(l^) above we should choose as fitting function 

G^{f, hi) = A(f)/if/'' + 5(f)/i,'/'' (59) 

(the notation f indicates that we are fitting at a fixed value of \rij\). 
However keeping only the first two terms in the scaling function of the 
correlator is a too rough approximation. Within the range of values of 
hi that we studied and with the precision that we obtained, this choice 
in general led to very high values. Fortunately, thanks to the IRS 
analysis we know the functional form of the next to leading orders in 
the hi expansion of the correlator. It turns out that adding one more 
term is enough for the correlators ranging from 1 to 7 lattice spacings 
both along the axis and the diagonal directions, while for the remaining 
correlators we have to go up to the fourth term in the expansion. The 
general form of the scaling function at this order turns out to be: 

G^{f, hi) = A(f)/if/'' + B{f)hf/'' + Cif)hf + D{f)hf (60) 
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(see for a discussion of this scaling function and the origin of the 
two new exponents hf^^^^ and h^'^^^^). In this way for all the values of r 
we found good confidence levels. The values of A{r) that we obtained 
in this way can be found in table |^. Notice that the data in each fit 
are completely uncorrelated since they belong to different simulations. 

3] As a last step we address the r dependence of the function A{r). Accord- 
ing to eq.(l^) we expect the following behaviour: 

A{r) , a^e) 

r 



Fixing ai to its known IRS value we end up with a one parameter fit 
which in both cases of = and = 7r/4 has a very good confidence 
level. Our final results are: 

• axis correlator 

02(0) = (-0.062 ±0.004) (62) 

• diagonal correlator 

a2(7r/4) = (-0.014 ±0.007) (63) 

These results, and in particular the fact that we have different corrections 
in the two directions, unambiguously show that the contribution is at least 
partially due to a spin 4 operator. The natural candidate for this role is 
again the ± operator. Our result also suggests that a scalar operator is 
present in the game. We see two possible reasons for this contribution. The 
first is that, since we are working in the IRS framework we are compelled 
to use the "continuum limit reference frame". As we have seen in sect. p.l.5| 
this choice induces a TT term and consequently a scalar type correction. 
The second is that one cannot exclude the presence in the magnetic scaling 
field of a (rotationally invariant) momentum dependent term, which would 
manifest itself exactly as a scalar correction of the type that we observe^]. 
This non trivial possibility has been recently discussed in (see the note 



at pag.8161) to explain the corrections of order (where t is the reduced 
temperature) in the second moment of the spin spin correlator. Most proba- 
bly both mechanisms are at work in the present case. We cannot distinguish 
between them within the precision of our data. 

^We thank A.Pelissetto for pointing to us this interesting possibihty. 
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Table 2: The lattice spacing N is intended along the axis for A"-^'^'^ and along 
the diagonal for A'^'"^^ . 



N 






1 


(0.581 ±0.003) 


(0.823 ±0.003) 


2 


(1.048 ±0.004) 


(1.396 ±0.004) 


3 


(1.446 ±0.005) 


(1.894 ±0.005) 


4 


(1.802 ±0.005) 


(2.354 ±0.007) 


5 


(2.134 ±0.006) 


(2.777 ±0.008) 


6 


(2.448 ± 0.006) 


(3.19 ±0.01) 


7 


(2.746 ± 0.007) 


(3.58 ±0.01) 


8 


(3.039 ±0.009) 


(3.96 ±0.01) 


9 


(3.316 ±0.009) 


(4.32 ±0.03) 


10 


(3.58 ±0.01) 


(4.67 ±0.04) 



6 Concluding remarks 

The role of the irrelevant operators in the two dimensional Ising model has 
attracted much interest in these last months due to the results on the mag- 



netic susceptibility at /i = recently reported in |T0|. While it is by now 
clear that contributions due to irrelevant operators are present in the free 
energy of the Ising model aih = Q nothing is known on the case in which the 
magnetic perturbation of the critical model is chosen. Moreover the charac- 
terization of these irrelevant operators and, possibly, their identification with 
quasiprimary fields of the Ising CFT is still an open problem. In this paper 
we tried to make some progress in this direction. We studied the corrections 
induced by the presence of irrelevant operators in the spin-spin correlator of 
the 2d Ising model in presence of an external magnetic field. 
Our main results are: 

• The 1/r^ corrections which are present in the correlator at the critical 
point survive unchanged also if the magnetic field is switched on. The h 
independent part of this contribution is much larger than our statistical 
errors and can be observed very precisely. It is completely due to the 
spin 4 irrelevant operator. 

• Several new terms appear when the magnetic field is switched on. By 
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using standard pcrturbative methods one can evaluate the amplitude of 
the corrections that they induce in the spin- spin correlator. In general 
these terms are very small. However in one case the amplitude of 
the expected correction turned out to be definitely larger than our 
statistical errors. It is the case of the term proportional to 

^ ■ (64) 



|^|l/4 \ J.2 j 

for which we could indeed observe and measure such correction with 
good confidence. The comparison between the two values of this cor- 
rection for the correlator along the diagonal and the one along the 
axis shows that, besides the expected contribution due to the spin- 
4 (T^ + T^) operator, there is also a scalar term which may have a 
twofold origin. It could be due the TT operator, but it could also be 
the signature of a momentum dependent term in the magnetic scaling 
field of the model. 

Other contributions due to these or other irrelevant operators are for the 
moment beyond our resolution, but it could be possible to observe them in 
the next generation of simulations. 
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A Determination of pi 

The simplest way to fix the value of pi is to look at the expectation value 
{(Ti). We expect: 

{a,)=RM+Pi{^) (A.l) 



from which we have (using the notations of ||27|| ) 

1 23 

(ai) = R^Auh-^ + PiAeH-^ (A.2) 



from 1^ we see that 



{ai) = ^4/.!^ + (A.3) 



Using the numerical values reported in i.e. 

A) = 0.9927995.... (A.4) 
4,2 ~ 0.021... (A.5) 

we obtain 

pi ~ 0.0345 (A.6) 



Acknowledgements We thank M. Hasenbusch, A. Pelissetto and E. Vi- 
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